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ABSTRACT 

We consider a model for gamma-ray bursts (GRBs) from high-velocity 
neutron stars in the galactic halo. In this model, bursters are born in the 
galactic disk with large recoil velocities Vr, and GRBs are beamed to within 
emission cones of half-angle 0^ centered on Vp. We describe scenarios for 
magnetically-channeled GRBs that have such beaming characteristics. We then 
make detailed comparisons of this halo beaming model (HBM) to data from 
the 3rd BATSE Catalog and from the Pioneer Venus Orbiter experiment, for 
both GRB intensity and angular position distributions. Acceptable fits to 
observations of over 1000 bursts are obtained for 06 = 15° — 30° and for a 
BATSE sampling depth of D ~ 180 kpc, which corresponds to a peak burst 
luminosity of ~ 10^'^ ergs s^^. Present data favor a truly isotropic (cosmological) 
model over the HBM, but not by a statistically compelling margin ( ^ 2o"). 

The HBM makes the distinctive prediction that the galactocentric quadrupole 
moment (cos^Q) — 1/3 for bright, nearby GRBs is large, even though the 
dipole moment (cos6) remains near zero. Bursters born in nearby external 
galaxies, such as M31, are almost entirely undetectable in the HBM because 
of misdirected beaming. We analyze several refinements of the basic HBM: 
gamma-ray intensities that vary with angle from the beam axis; non-standard- 
candle GRB luminosity functions; and models including a subset of bursters 
that do not escape from the galaxy. We also discuss the energy budgets for 
the bursters, the origins of their recoils, and the physics of burst beaming and 
alignment. One possible physical model is based on the magnetar model of 
soft gamma repeaters (SGRs). Empirical bounds on the rate of formation and 
peculiar velocities of SGRs imply that there exist ~ 10^ to ~ 10^ aged SGRs in 
the galactic halo within a distance of 100 kpc. 
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The HBM gives an acceptable fit to observations only if it satisfies some 
special conditions (0^ ~ 20°, uniform bursting rate) which are possible, but 
for which there are no clear and compelling theoretical justifications. The 
cosmological burst hypothesis is more generic and thus more attractive in this 
sense. 

Subject headings: gamma rays: bursts — stars: magnetic fields — stars: neutron 
— galaxy: halo 
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1. Introduction 

The Burst and Transient Source Experiment (BATSE) on the NASA Compton Gamma- 
Ray Observatory has revealed a nearly isotropic but inhomogeneous sky distribution of 
gamma-ray bursts (Meegan et al. 1992; pVIeegan et al. 19961) . An additional major constraint 
on GRB theories, revealed by the Pioneer Venus Orbiter (PVO) experiment, is that the 
intensity distribution of bright bursts is consistent with a locally-uniform density of sources 
in Euclidean space ( [benimore et al. 1993| ). These facts are simply and generically accounted 
for if the faintest observed GRBs come from cosmological distances (e.g., Paczyhski 1995). 
An alternative possibility, that we will focus on here, is that the bursts come from the 
extended halo or corona of our galaxy. Galactic halo models for GRBs were first considered 
and tested against data by Fishman et al. (1978) and Jennings & White (1980). Reasons 
for favoring halo models have been summarized by Lamb (1995). 

In galactic halo GRB models, it is often supposed that the bursters are neutron stars, 
since such compact stars are conjectured to be capable of producing intense fluxes of hard, 
non-thermal photons; while their small size could drive burst variability on submillisecond 
time scales. Other reasons for favoring neutron stars include observations of spectral lines, 
which are at present controversial, and a possible connection of classic GRBs with the 
March 5, 1979 gamma ray burster (e.g., puncan fc 'I'hompson 19W2 , hereafter DT92) which 



was localized to an angular position lying within a young supernova remnant ( |Gline et 



al. 1982| ; [Rothschild, Kulkarni fc Lingenfelter 1994]) . The displacement of this burster from 
the center of the supernova remnant indicates that it acquired a velocity ~ 1000 km s~^ at 
birth (DT92). Neutron stars with such velocities will escape the galactic disk, and move 
into the halo. Recent analysis of the 1979 March 5th event lends support to its association 
with classical GRBs ( |Fenimore, Klebesadel, fc Laros 1996 , hereafter FKL). In addition. 



many high- velocity radio pulsars {V > 500 km s ^) have been observed ( [Lyne fc Lorimer 
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1994; Ibrail 1996|) . 



For these reasons, we will focus on theories of GRBs from high-velocity neutron 
stars (HVNSs) born in the galactic disk, as first suggested by Shklovskii & Mitrofanov 
(1985). In particular, we will consider the "halo beaming model" (HBM) proposed earlier 
( Puncan, Li, fc Thompson 1993| , hereafter DLT; |Li, Duncan, fc Thompson 1994| , LDT; Li 



& Duncan, 1996a, 1996b; Bulik & Lamb 1996). Alternative models for GRBs from HVNSs 
in the galactic halo invoke a delayed onset of bursting activity at a time ~ 10^ years after 
birth in the disk ( [Li fc Dermer 19"9^ , hereafter LD92), which has been applied in several 
different physical contexts ( [Colgate fc Leonard 199^ ; [Lamb, Bulik fc (Joppi 19^B| |Woosley| 
|fc Herant 1996| ); and weakly-bound bursters orbiting in a nonspherical galactic potential 
( Podsiadlowski, Rees fc Ruderman 1995| , hereafter PRR). Models for bursters born in the 



galactic halo ( pichler fc Silk 1992| ; [Hartmann 1992| ; ^alpeter fc Wasserman 1993| ; [Wassernian 



fc Salpeter 1994|) , or the Magellanic Clouds ( [Fabian fc Posiadlowski 1993| ), are beyond the 



scope of this paper. 

In §0, we will describe our basic model for the galactic population of bursters, and 
briefly review its physical motivations. In §^ we present our basic Monte Carlo model 
results. We discuss how these results arise, and what they could imply for halo GRB 
theories. We also explore the model's sensitivity to different beaming angles. In §^ we make 
detailed quantitative comparisons of the model with BATSE and PVO data. In particular, 
we study moments of the angular position distribution in bright subsets of observed GRBs, 
which is potentially a sensitive model discriminant. In §|^ we estimate the GRB repetition 
rate and the energy requirements for bursters. We also discuss two candidate power sources 
that could satisfy these requirements: magnetic energy and accretion energy. Several 
refinements of the basic HBM are investigated in §^, namely: gamma-ray intensities that 
vary with angular position within the beam, GRB luminosity functions, and models with 
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a significant subset of bursters on bound orbits in tlie galactic lialo. In §0 we give our 
conclusions and outline future observational tests. 

Note that we consider possible physical mechanisms for halo GRBs, involving neutron 
stars with unusually strong magnetic fields, in §2.3, §2.4, §5, §6.3, and Appendix B. These 
sections could be read separately from the rest of the paper, since they might apply in a 
non-beamed model context (cf. §7.2). 

2. The Halo Beaming Model: Physical Motivations 

2.1. Model Assumptions 

We will assume that GRBs are emitted by HVNSs, which emanate from the galactic 
disk like a "wind" extending into the galactic corona (|Shklovskii fc Mitrofanov 1985|) . 



The HBM does not require that the bursters or their peculiar velocities have any favored 
orientation in a galactic coordinate frame, but it does invoke the physically-plausible 
condition that gamma-rays are produced only within a cone of angular radius (pb about the 
star's magnetic axis ±7?. The magnetic axis is furthermore assumed to be roughly aligned 
(within ~ 20°) with the stellar recoil velocity Vr. We discuss the physics of such beaming 
and alignment in §^72. 



The particular version of HBM which we will analyze quantitatively in §^ and §§ below 
has the following simple properties: [1] bursters are born at positions distributed like young 
Pop. I stars in the galactic disk; [2] with randomly directed recoils Vr = 1000 km s^^; [3] 
they emit GRBs at a constant rate, with [4] constant luminosity, and [5] the gamma ray 
emission is beamed parallel and anti-parallel to Vr, within an angular radius 0b that we 
will vary. 

GRB beaming in a direction correlated with Vr makes the observable burst distribution 
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comply with the BATSE dipole isotropy and the PVO brightness-distribution constraints, 
for the following reasons (see Figure ||). Since the HVNSs are freely-streaming out of the 
galaxy, their mean density diminishes with distance from the galactic center as n ~ r~^. 
However, the fraction of escaping bursters that are potentially detectable at Earth increases 
in proportion to the transverse area of their beaming cones, ~ r^. These two trends cancell, 
making the effective number of bursters increase linearly with the sampled volume (i.e., an 
apparent "constant density" of detectable stars) within a "core radius" Re ~ Ro/(pb, where 
Ro ~ 8.5 kpc is a galactic disk dimension. In the HBM, this produces the the observed 
"homogeneous" distribution of bright GRBs found by the PVO experiment. Furthermore, 
since most bursters are undetectable when they are at small r, the observable dipole 
anisotropy of GRB positions in the direction of the galactic center is greatly reduced. At 
distances larger than Re, all bursters are detectable at Earth (or nearly all; see exceptions 
discussed in §3), and the n oc r~^ free-streaming fall-off of burster density prevails, 
accounting for the "boundedness" {{V/Vmax) < 0.5) found by BATSE. 

A more detailed illustration of the geometrical effect of burst beaming is given in 
Figure ^. To understand this figure, consider for a moment an idealized model in which all 
stars are born precisely at the galactic center (GC), and move out in random directions on 
straight-line trajectories, each emitting GRBs into a cone of half-opening angle (ph around 
its velocity vector. Bursters lying outside the circles in Figure |^, or within the lens-shaped 
intersection of circles between the GC and Earth, are then the only ones that can be 
detected at Earth. There is a large "zone of avoidance" (ZOA), within which all bursters 
are invisible. Figure ^ actually shows only the 2-D cross-section of this ZOA. The true 
ZOA is the volume of revolution of the pictured shape about the line between Earth and 
the galactic center. The boundary in Figure || at which stars just become observable (edge 
of the ZOA) is part of a circle with radius Z^sun/sin0fe, where Dg^^ = 8.5 kpc is the distance 
from Sun to the galactic center. 



In the realistic situation, bursters are born throughout the galactic disk, with a peak 
birth rate at ~ 4-5 kpc from the center, where the greatest concentration of Population I 
stars are located ([van der Kruit 1987|) . Each birthplace in the disk then has its own ZOA, 
scaling up or down in linear size with the distance between the birthplace and Earth. One 
must add the weighted distributions of detectable bursts together to get the total observable 
GRB distribution. There is no simple, analytic way to do this, so in what follows we will use 
Monte Carlo methods. We will also calculate realistic trajectories in the galactic potential, 
rather than assuming straight lines. 

2.2. Physics of Beaming and Alignment 

Gamma-ray bursts from high-S environments tend to be beamed along field lines 
because of the transverse pair-creation opacity (e.g., [Riffert, Meszaros fc Bagoly 1989 ; 



Ho, Epstein fc Fenimore 19901) and because gamma rays produced by such mechanisms 



as curvature radiation and Compton upscattering are strongly beamed along field lines 



(e.g., [Sturrock 1986| Dermer 1990 and references therein). Even locally-isotropic emissions 
are strongly beamed when they occur in a relativistic outflow channeled by a magnetic 
fleld (e.g., Yi 1993). Gamma emissions induced by sheared Alfven waves in a neutron 
star magnetosphere can also be highly beamed (Melia & Fatuzzo 1991; Fatuzzo & Melia 
1993). Beaming obviates the 7-7 opacity limit for GRBs (Krolik & Pier 1991; Baring 1993; 
Fenimore Epstein & Ho 1993). 

The HBM invokes the additional (assumed) property that burster recoil velocities are 
approximately aligned with magnetic dipole — and hence burst beaming — axi: Vr \\ Jt. 
This requires that the rotation axis Q is also roughly aligned with Jt, at least to within the 
burst opening angle (pb. 
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Is Vr II 7^ a plausible assumption? Any recoil mechanism which imparts impulses to 
the stellar surface with a coherence time longer than the rotation period of the star will 
tend to make Vr \\ ±^ because of "rotational averaging." One such recoil mechanism is 
anisotropic neutrino emission during the first ~ 10 s after formation (§2.4 below). If it 
is also true that 7^ || ±i7 (at least to within the beaming angle 0b) then the alignment 
condition of the HBM is satisfied. 

Such near-alignment between jt and ± fi is expected if the neutron star magnetic field 
is generated by large-scale dynamo action, as in familiar stellar and planetary dynamos. 
The "magnetar" model, described below (§2.3), is one scenario involving such a dynamo. 
Vacuum magnetic torques increase alignment ( [Michel fc Goldwire 197(J| ; Davis fc Goldstein 



1970| ) after the star spins down past the "death line" at which magnetospheric currents 
are quenched, at age t^, ~ 10^ (-Bdjpo^e/lO^^ G)^ yrs (Chen & Ruderman 1993). Before this 
point is reached, currents might exert counter-aligning torques of comparable strength (e.g., 
McKinnon 1993 and references therein). 

Note that the GRB beaming angle (pb is the sum of the intrinsic beaming angle due 
to gamma-ray opacity and radiative effects, and the r.m.s. scatter in the angle between 
Q and Jt. Gamma production mechanisms such as curvature radiation and Compton 
upscattering could operate at significant altitudes in the magnetosphere, where field hues 
are tilted with respect to 'Jt, yielding larger values of (pb than one would estimate under the 
assumption of near-polar gamma emission (as adopted by [Rihert, IVleszaros fc Bagoly 1989 



Ho, Epstein fc I'enimore 1990| ). Values as large as (pb ~ 20° are possible. 
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2.3. Magnetically-powered Bursts? 

The HBM might apply in a variety of different physical contexts. Here we will briefly 
discuss the "magnetar" model for GRBs (DT92). The magnetar idea has previously been 
invoked in models of soft gamma repeater bursts ([Thompson fc Duncan 1995| , hereafter 
TD95) as well as in models of continuous X-ray emissions from SGRs and from anomalous 
X-ray pulsars like 1E2259+586 ([Thompson fc Duncan 1996| , TD96). 



Magnetars are hypothetical neutron stars which are born with dipole fields in excess 
of Bq = me^c" /eh = 4.4 x 10^^ G. How might such stars form? Neutron stars are hot 
and convective during the first ~ 10 seconds after they form. They also undergo strong 
differential rotation (TD93). Large-scale dynamo action might operate efficiently during 
this time for neutron stars with mean rotation periods below a threshold comparable to the 
convective overturn time, as predicted by a-Q dynamo theory. There is evidence for such 
a threshold in magnetically-active main-sequence stars (e.g., Simon 1990). This threshold 
effect could give rise to a bimodal population of neutron stars ("pulsars and magnetars") 
with a factor ^ 10^ difference in field strength (DT92). A more detailed explanation of 
this magnetar formation hypothesis is given by Duncan & Thompson (1996, DT96). Seven 
distinct estimates of the field in the bursting neutron star SGR 0526—66 seem to indicate 
B > 10^^ G (TD95).g 



■^Katz (1996) has raised some interesting questions about theories of SGR burst emissions. 
Our views on these points are given in TD95. In particular, equation (2) of TD95 explains 
why we believe that plasma-confining magnetic fields ^ 10^^ G are necessary in SGR 0526- 
66; §3.2 of TD95 explains why we do not favor a model for SGR bursts involving a pure pair 
plasma confined by a ~ 10^^ G field, as favored by Katz; and §7.3.2 explains why we do not 
favor a hydrostatic model in which the magnetic Eddington limit operates without confining 
magnetic stresses, which Katz incorrectly attributes to us. Another good argument against 
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Magnetars spin down too rapidly to be observed as radio pulsars. DT92 conjectured 
that GRBs are flare-like reconnection events in a galactic halo population of such 
strongly-magnetized neutron stars. We will show in §4.5 that the available magnetic energy 
is sufficient to power observed GRBs, in the context of the HEM (see also PRR). 

Magnetic reconnection is a theoretically advantageous energy source for GRBs, for 
several reasons. Gravitational and nuclear energy releases generally occur in bulk baryonic 
matter, which has many degrees of freedom into which energy can thermalized and degraded 
via adiabatic expansion (e.g., |Piran fc Shemi 1993| ). Flares in a neutron star magnetosphere. 



on the other hand, can be "clean," exciting only photon and pair degrees of freedom to a 
first approximation. Large scale electromotive forces induced by reconnection will generally 
accelerate pairs, and may produce hard, non-thermal spectra via Compton upscattering, 
synchrotron emission, and curvature radiation ( ^turrock 1986|) . Indeed, GRBs have many 
qualitative similarities to stellar flares. Both are transient energy releases, with chaotic 
variability over a wide range of time scales; and both have spectrally-hard non-thermal 
components (e.g.. Murphy et al. 1993; Ramaty & Mandzhavidze 1993). The similarities are 
so strong that a minor fraction of bursts in the BATSE catalog might actually be intense 
events from nearby, non- neutron, flare stars (Liang & Li 1993). Stochastic avalanche models 
give a good fit to GRB time profiles (Stern & Svensson 1996) and to many properties of 
solar flares (Lu et al. 1993). 

If the Vr-Bdipoie correlation observed in some samples of radiopulsars (e.g., Cordes 
1986; StoUman & van den Heuvel 1986) is applicable for all neutron stars, then the mean 
magnetar recoil velocity would be ^10 times larger than the mean for radiopulsars (but 
see Lorimer, Lyne and Anderson 1995). Several mechanisms predict unusually large recoils 
for magnetars, sufficient to propel them into the galactic halo, as explained in §2.4 below. 



hydrostatic models, based on some novel physics, is given by Miller (1995). 
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Because of the dynamo origins of magnetar fields, one expects that initially fl , ± K ^ 
and ±7? are roughly aligned in niagnetars.0 As the star spins down, niagnetospheric 
currents might drive some degree of misalignment; however, once the star spins down past 
the "death line" at time ~ 10^ {Bdipoie/ 10^^ G)~^ yrs ( |Chen fc Ruderman 1993]) , further 



spindown certainly enforces alignment ([Michel fc Goldwire 1970] ; Pavis fc Goldstein 19701) 



since the magnetostatic stellar distortion in magnetars is large enough to damp nutations of 
Q about jt ( |Goldreich 1970| , DLT). Incidentally, this is probably not true in old, spun-down 
{Prot > 4 s) radiopulsars. We conclude that Vr is likely to be roughly aligned with ±jt in 
old magnetars, where ±~jt is also the axis of beamed gamma emissions (DLT).0 

This scenario for classic GRBs assumes that magnetars retain strong dipole magnetic 
fields. Large-scale magnetic instabilities that reduce Bdipoie (Flowers & Ruderman 1979) 
might be suppressed by toroidal field components in the stellar interior, as expected for 
fields generated via a-Vt dynamos, or perhaps by other mechanisms (TD93 §14.2). Note 
that a dipole field anchored in the stably-stratified liquid interior of a magnetar cannot 
be greatly distorted by spindown-induced crustal tectonic drift (Ruderman 1991) because 
magnetic stresses dominate tensile stresses in the crust. This might differ markedly from 
the circumstance in ordinary radiopulsars (Ruderman 1991). The interactions of vortex 



^Magnetically-induced neutrino starspots, as discussed in §2.4, will have coherence 
times Tapot ^ Prot when the convective overturn time satisfies Tcon ~ Prot, thus "rotational 
averaging" implies that Vr \\ ±^ (where Vl = 2n/Prot)', and Vl \\ ±~jt because the dipole 
field is generated by an a-Q dynamo. 

^Fenimore, Klebesadel & Laros (1996) have suggested that an initial period of (spindown- 
induced) misalignment occurs during the SGR phase, allowing bursts hke the 1979 March 
5th event to be observed. Subsequent alignment by vacuum magnetic torques occurs when 
classic GRB emissions begin. 
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lines in the neutron superfluid with superconducting flux tubes are also probably irrelevant 
for young magnetars, because the interior field is probably strong enough to suppress 
superconductivity. A young magnetar's field evolves predominantly via ambipolar diffusion 
in the interior and Hall fracturing in the crust (TD96). As the interior field decays, 
superconductivity will appear in the mantle and perhaps the core; but this probably 
happens only after rapid spindown has driven most superfluid vortex lines out of the 
interior. Further discussion of magnetar evolution is given by PRR. 

2.4. Neutrino Magnetic Recoils 

A dipole anisotropy of only ~ 0.03 in the neutrino emission from a young, hot neutron 
star would impart a recoil velocity ~ 1000 km s~^ to the star (Chugai 1984). Here we 
briefly review mechanisms whereby neutrino anisotropics could be magnetically induced. 
Purely hydrodynamic (non-magnetic) mechanisms for producing neutron star recoils have 
been proposed by Janka & Miiller (1994), Shimizu, Yamada & Sato (1994), Burrows & 
Hayes (1996), and in references quoted therein. 

A strong magnetic field affects neutrino emissions from the beta processes n -^ p^ e~ u 
and p"*" e~ — *> n z/, as calculated by Dorofeev, Rodionov & Ternov (1985), and from 
neutrino scattering processes as calculated by Vilenkin (1995). Dorofeev et al. and Vilenkin 
furthermore estimated the neutron star recoils resulting from these processes, in the 
uniform-field idealization. This is a macroscopic manifestion of parity-nonconservation in 
the weak interactions.^ 



^ Another neutrino magnetic recoil estimate, by Bisnovatyi-Kogan (1996), makes several 
assumptions: (1) the neutron star magnetic field undergoes simple linear winding driven 
only by differential rotation (i.e., ignoring all mixing motions); (2) the winded field is 
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In a realistically non-uniform magnetic field, the back-reaction of magnetic stresses on 
convective energy flow, along with the neutrino opacity variations outlined above, would 
give rise to "neutrino starspots," analogous to sunspots, and thus produce neutron star 
recoils of a magnitude estimated in DT92 and §13 of TD93. 

These calculations are based on standard Weinberg-Salam weak interaction theory. 
If neutrinos have mass, on the other hand, a strong magnetic field will shift resonant 
flavor- changing (MSW) oscillations, thereby also inducing an anisotropy in the emergent 
neutrino flux from a nascent neutron star (Kusenko & Segre 1996). 

All of these mechanisms produce recoils which vary directly — usually linearly — with 
the mean magnetic field strength: Vr cc B. Thus if any of these mechanisms dominate, 
one would expect much larger recoils for magnetars than for pulsars. Several non-neutrino 
recoil mechanisms which also might operate more efficiently in magnetars than in pulsars 
were discussed by DT92.[] 

To estimate neutrino magnetic recoils realistically, one must know the strength. 



largely cancelled by a postulated pre-existing toroidal field in one hemisphere only; (3) 
the free neutron decay rate in a magnetic field applies. These assumptions are questionable. 
In particular, realistic estimates of the beta neutrino emissivity within a nascent neutron 
star should take into account the presence of magnetized, degenerate, relativistic electrons 
(Dorofeev, Rodionov & Ternov 1985). Dorofeev et al.'s work was itself an improvement on 
Chugai's pioneering calculation, which did not treat the energy levels of the electrons in a 
B ^ Bq field correctly. 

^A star with an oblique, off-center dipole field experiences a spindown- driven recoil 
(Harrison & Tadmaru 1975), but this is unlikely to produce velocities as large as ~ 1000 km 
s~^, as explained in DT92. 
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coherence length and coherence time of the magnetic field during the epoch of maximum 
neutrino luminosity. This is not possible at present. However, if the magnetic field 
approaches equipartition with the free energy of differential rotation and/or the convective 
fiuid mixing — which is strongly in the MHD limit — then i? x^ 10^^ G at the neutrinosphere 
(TD93). Such strong fields could be present when most of the neutrino energy is radiated 
away, even if the surface dipole field, which is frozen-in tens of seconds later by the onset 
of stable stratification in the liquid interior (Goldreich & Reisenegger 1992), is smaller by a 
factor of ~ 10 or ~ 30. Recoils Vr ^ 10^ km s^^ are plausible. 

3. Galactic Halo Model Results 

By integrating a large number (~ 10®) of neutron star trajectories in a realistic galactic 
potential, we have derived the sky distribution of HVNSs. In the model described here (§3 
and §4), all stars have |Vr| = 1000 km s~^; thus they move in nearly straight lines and 
eventually escape from the galaxy (but see §6.3). 

We used a Monte Carlo code developed by Li & Dermer (1992), which in turn follows 
many of the prescriptions of Paczyhski (1990) and Hartmann, Epstein & Woosley (1990). 
In particular, we adopt van der Kruit's (1987) model for the spatial distribution of young 
Pop. I stars which spawn neutron stars in the galactic disk. The trajectories of ~ 10® 
neutron stars were numerically integrated in a realistic galactic potential (Miyamoto & 
Nagai 1975) including a dark halo cutoff at radius Rh = 70 kpc. The escape velocity from 



'^Lamb, Bulik, & Coppi (1996) have argued that the galactic potential of Miyamoto 
& Nagai (1975) implies an unrealistically extended disk which tends to 'focus' neutron star 
orbits into the disk plane. We find, however, that this focusing effect is negligible for randomly 
oriented initial recoil velocities of magnitude Vr ^ 1000 km s~^. This is because the disk is 
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the center of this model potential is Vesc ~ 600 km s~^. This increases only logarithmically 
with Rh, reaching 800 km s~^ for Rh = 200 kpc, thus our results are probably insensitive 
to Rh over its plausible range (LD92). Deviations from sphericity in the dark halo were 
neglected (but see PRR). We also do not include the potential of M31 since stars with 
Vr ~ 10^ km s'^ are negligibly perturbed by M31 within the BATSE sampling depth we 
found of < 200 kpc (§4.1). 

In Figure § we show HBM numerical results for several angular statistics, plotted as 
functions of sampling depth D about the Earth. That is, the figure shows cumulative 
angular statistics for all detectable bursters at distances from Earth that are less than or 
equal to the value of D on the horizontal axis. The topmost plot of Figure ^ shows the 
galactocentric dipole moment (cosB) , where G is the angle between a burst and the galactic 
center. The second plot shows the disk-like quadrupole, (sin^6) — 1/3, where b is galactic 
latitude; (sin^6) < 1/3 implies that sources are concentrated toward the disk. The third 
plot shows the galactocentric quadrupole, (cos^G) — 1/3. The bottom plot of Figure |] 
shows {V/Vmax), a statistic which is related to the slope of the cumulative log Af — log P 
brightness distribution, as explained below. 

Within each subplot of Figure ^ the various lines correspond to different beaming 
angles (pb as described in the figure caption. We have cut off the (pb = 10° and 5° curves for 
D < 10 kpc because our Monte Carlo sampling statistics at smaller D are too poor. For 
example, the number of detectable stars aX D < 3 kpc in our Monte Carlo model is only 
~ 55 for (pb = 5°, whereas it is ~ 10^ for in the unbeamed case {(pb = 90°). 

Increasing the sampling depth D (moving to the right in Figure |^) is tantamount to 
including fainter and fainter bursts. Eventually the BATSE sampling depth is reached; at 



only a very minor part of the total mass of the Galaxy. 
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this point, if the model is to fit observations, all the plotted statistics must simultaneously 
match BATSE values to within observational uncertainty. The most recent published 
BATSE results^ are plotted in Figure |] at a value D = 180 kpc. Before discussing these 
results, we must explain the significance of {V/Vmax)- 

The V/Vjnax statistic was invented for the study of the quasars (Schmidt 1968), and 
was first applied to GRBs by Schmidt, Higdon & Hueter (1988). For scintillation counter 
experiments, {V/Vmax) is the average of (C/Cmin)^^^'^ over all bursts, where C is the 
peak counts (in a set time interval) and Cmin is the threshold for detection. Cmin can 
vary with background noise and other effects (e.g., "overwrites"); however, for a source 
population that is uniformly distributed in static Euclidean space, {V/Vmax) is equal to 0.5 
regardless of how the threshold varies; and (V/Vmax) < 0.5 indicates that the density of 
bursters diminishes with distance, for standard-candle sources. Because Cmin is variable, 
the (V/Vmax) statistic is useful in the study of inhomogeneous data sets only when trying 
to answer the yes/no question: "Is the observed brightness distribution consistent with 
a uniform density of sources distributed in Euclidean space?" (e.g.. Band 1992; Petrosian 
1993). 

In Figure ^, we use (V/Vmax) in a purely illustrative way. Our model values of (V/Vmax) 
are ideal values that would be found by an instrument with a uniform detection threshold, 
i.e., no variations in the noise or the threshold settings. Since the Cmin values in the BATSE 
catalog are not highly variable (Meegan et al. 1996) preliminary comparisons with BATSE, 



^ Results from the third BATSE catalog for 1122 bursts (Meegan et al. 1996) are as follows: 
{cosQ)data = —0.002 ± 0.017. There is a concurrent anisotropy in the BATSE sky coverage 
of {cosQ)sky = —0.013. Thus the best measure of the burst population dipole moment is 
(cosB) = {cose)data " (cos0),fcj, = +0.011 ± 0.017. Similarly, (sin^ft) - 1/3 = 0.002 ± 0.009, 
and for 657 bursts, (V/Vmax) = 0.33 ± 0.01. 
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as in the bottom panel of Figure |^, will not lead us astray. However, the most accurate 
way to compare burst observations with theory is to fit the observed distribution of peak 
photon fluxes to Monte Carlo models which have been realistically filtered for detection 
incompleteness (e.g., Lubin & Wijers 1993). This is what we do in our actual statistical 
comparisons of the HBM with the BATSE catalog (§4). 

Note that the model curves for angular statistics in Figure |^ implicitly assume a 
detector with uniform sky coverage; i.e., we have idealized that the detector is equally 
capable of detecting bursts from any location on the celestial sphere. We have corrected for 
the imperfect sky coverage of BATSE by shifting the data points in Figure ^ appropriately 
(see footnote 9). In §^we will take the opposite approach: using the raw BATSE data and 
filtering the Monte Carlo model to take into account BATSE's imperfect sky coverage. 

What can be learned from Figure ^ Beaming evidently increases (V/Vmax) to nearly 
0.5 for nearby bursters, D ~ 30 kpc. This can make the model satisfy PVO constraints, as 
we show quantitatively below. Beaming also evidently reduces (cosG). The 0f, = 20° case 
fits BATSE observations over an appreciable range of sampling depths D > 100 kpc. 

Note that the largest deviations from isotropy in Figure ^ occur for bright (i.e., 
small--D) subsets of the observable bursts. A distinctive signature of the HBM is that 
(cos^0) — 1/3 is positive for bright bursts. This can be understood as follows. The (bright) 
bursters which become visible at Earth before reaching the remote galactic halo are the ones 
which happen to be born with beaming axis 'Jt (and hence, with Vr) pointed approximately 
toward or away from Earth (to within ~ 0^). Because most bursters are born within the 
Solar circle in the galactic disk, the brightest ones tend to be seen in that direction and 
toward the galactic anticenter, if they have already moved past the Earth on their way 
out of the galaxy. Thus (cos^O) > 1/3 for bright bursts, while the dipole moment (cos0) 
remains small. 
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AU burst sources can be seen at distances D ^ Ro/<Pb in the (ph = {90°, 40° and 
20°} cases of Figure |^, thus in a first approximation, one would expect (sin^6) to reach 
the isotropic value of 1/3 asymptotically at large D. In fact, a mild degree of disklike 
quadrupole anisotropy, (sin^6) = 0.31, remains in these models at large D. This is due to 
the galactic disk rotation (with velocity V^ ~ 220 km s~^), which makes bursters tend to 
stream out at lower galactic latitudes, diminishing the density of bursters near the galactic 
poles compared to the equator at a given distance r from the galactic center (LD92). This 
effect is not dependent on beaming. It depends only on the ratio Vq/Vt, with (sin^6) -^ 1/3 
at large D as V^/K —>■ 0. 

Galactic rotation also causes the total initial velocity vector Tv + Vq to be inclined by 
as much as as 

C-.a.cta„(VV,)«12"(^555^)" (1) 

with respect to the axis of GRB emission, Tv/K- Thus if (pb < 0™"^, some stars never 
become visible from Earth. Since this happens preferentially at high galactic latitudes, 
(sin^fo) is always substantially less than 1/3 in models with (ph < 0^"^. For example, in the 
case (pb = 10°, (sin^6) never exceeds 0.25. Only when (pb > 0™"^ are all stars visible at large 
D, causing (sin^6) to asymptote to the unbeamed, large-Z) value of 0.31. 

When (pb < 0™"^ most stars eventually become undetectable at Earth, 
but only after passing through a phase when they can be seen, namely, when 
D < 8.5kpc/[tan (0™"^' — 0^)] ^ 150 kpc. This is evident in the quadrupole moment plots 
of Figure |^. 

This leads us to consider bounds on the range of acceptable beaming angles (pb- 
The effective "core radius" in the HBM, is Re ^ Ro/(pb ^ 24(06/20°)"^ kpc. This value 
is confirmed by our numerical results if we define Re as the sampling depth at which 
(V/Vmax) = 0.46 close to the value found by PVO (Fenimore et al. 1993). Thus by equation 



-1 
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(0), 

^c-40V3(^j kpc, (2) 

where V3 = K/IO^ km s"^ Since BATSE bounds on (sin^6) imply (06/0^"^) > 1, the 
HBM cannot have a core radius larger than ~ 40 V3 kpc. The minimum core radius that 
is likely to fit the data is Re ~ 16 kpc (Paczyhski 1991), thus the range of allowed (pb is 
roughly 

12° V^^ <(l)b< 30° . (3) 

This agrees with Figure 0, where (pb = 20° give the best fit to BATSE data. 

Although Fig. 3 is for Vr = 1000 km s^^, note that better fits to the BATSE data can 
be obtained when Vr > 1000 km s^^. This reduces the effect of galactic rotation on the 
burster distribution, causing (sin^fe) to converge toward 1/3 at the faint end, and it also 
widens the acceptable range of (pb (eq. [Q). Some soft gamma repeaters have recoils as large 
as Vr ~ 2000 km s"^ (Appendix A). 

4. Model Comparisons with the 3B Catalog and PVO 



We now make detailed comparisons of the HBM with the 3rd BATSE catalog ( [Meegan 



et al. 199^ ), which we will refer to simply as "3B". Because the largest deviations of angular 



moments from isotropy in the HBM occur within bright subsets of bursts (small D in Figure 
^, we will explicitly consider the cumulative angular statistics as a function of observed 
peak photon fluxes, P. This is the basic approach that we previously applied to the flrst 
and second BATSE catalogs (LDT; LD96a). A preliminary discussion of 3B comparisons is 
given by LD96b; see also Bulik & Lamb (1996) for a complimentary statistical analysis. 
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4.1. The Log AT-Log P Distribution 

There are several issues that must be faced when doing statistical analyses of the 3B 
catalog. There are data gaps in in many bursts, and for many bursts the associated trigger 
efficiency (the probability that BATSE will detect a burst, as a function of peak flux) is less 
than one. Under these circumstances, we use the following criteria for selecting bursts from 
BATSE and PVO: (1) Overwriters in 3B are not included. (2) Only bursts from 3B with 
-^256 ^ 1 photons cm^^ s^^, whose trigger efficiency presumably is > 1, are used. Similarly, 
bursts from PVO with P250 ^ 25 photons cm"^ s^^ are used. (3) Only bursts for which 
(C'max/C*min)256 > 1 or (Cmax/C'min)io24 > 1, or both, are used. This applies to both 3B and 
PVO. 

This counting criterion reduces the number of BATSE bursts used in our study of the 
log AT-log P distribution to 241 out of the total number 1122 in 3B.Q This, however, also 
allows us to construct a relatively homogeneous merged catalog with the PVO data, since 
the PVO bursts are selected by very similar criteria (Fenimore et al. 1993). Thus we add 
181 PVO bursts to our analysis of the log A/" — log P distribution, making a total of 422. 
More importantly, by adding PVO we have increased the range of peak fluxes which we can 
flt by a factor of about 20. 

We have also done analyses by selecting bursts according to criteria (1) and (2), and 
ignoring criterion (3). This is somewhat justifled since the BATSE bursts for which P256 
and Cmax/Cmin are available all have (Cinax/Cmin)256 > 1 or (Cmax/Cmin)io24 > 1 when 
-P256 ^ 1 photons cm^^ s^^. This gives us a larger sample, with 437 bursts from BATSE 



^'^ We will use many more 3B bursts in our study of the GRB angular distribution, 
as described below. Note that a signiflcant fraction of 3B bursts do not have Cmax/C'min 
available, so we had to drop these bursts from the start. 
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and 198 from PVO, thus making a total of 635. The fitting parameters we obtained from 
this sample, especially the peak luminosity, are very similar to the values we obtained from 
fitting the 422 sample, however. Thus we will only present the results from the sample of 
422. 

The combined PVO and BATSE intensity distribution is shown in Figure ^. The HBM 
result for 0{, = 20° and Vr = 1000 km s~^ is also shown. We have used two free parameters 
in fitting these data to the HBM: the absolute number normalization (i.e., the number 
of GRBs per year) and the intrinsic peak photon emission rate A^. We find a best-fit 
"standard candle" value of 

N ~ 9.45 X 10^^ f ^-cos(pA pi^Q^^^g g-i (4) 

This fit has a minimum x^ — 6.9 with 12 degrees of freedom. The l-a confidence region of 
N, estimated as the range over which x^ increases by 2.3, is 8.81 — 11.05 x 10^^ photons 
s~^. The corresponding peak luminosity is 

i^ - 1.5 X 10^° (^^) ergs s-^ (5) 

if the average photon energy is ~ 100 keV. The corresponding BATSE sampling depth is 
-Dbat ~ 182 kpc for a limiting flux of ~ 0.3 photons cm^^ s^^. 

We have also fit the model to 3B alone, ignoring PVO constraints. This is shown in 
Figure ^. The best-fit peak luminosity turns out to be 

iV = 6.9x10^^ /^^cos0fe\ pj^Q^Qj^g g-i^ BATSE alone (6) 

with a 1-0" confidence region ranging from 4.7— 10.4 x 10^^ photons s~^. This is smaller than 
the number we got using the combined PVO and BATSE data, although the confidence 
regions do overlap. 
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4.2. The Angular Distributions 

In order to compare the HBM with the observed burst angular distributions, we have 
filtered our model results through BATSE's sky exposure table. All bursts from 3B with 
tabulated values of P256 are used in this comparison except overwriters, giving us 814 bursts 
in total. 

In Figure P, we plot measures of angular (an)isotropy in galactic-based coordinates for 
bursts from the 3B catalog in which the peak photon flux (P256) is greater than or equal to 
a given value, shown on the horizontal axis. Here (pb = 20°, and we have used the best-fit 
photon luminosity N of equation |[ 

The HBM curves (solid lines) make the distinctive prediction that (cos^ O) > 1/3, as 
explained in §^. Dashed lines show model predictions if gamma-ray bursts are isotropic on 
the sky at all intensity levels, as in the cosmological burster hypothesis. Note that both 
model curves (HBM and isotropy) are corrected for the incomplete BATSE sky coverage. 
The numerical (Monte Carlo) uncertainties on the HBM curves are very small on the scale 
of this graph: comparable to the width of the thick solid lines. 

Filled dots with error bars in Figure ^ are BATSE cumulative data points with 1-cr 
statistical uncertainties evaluated according to 



a 



((X)) = N-^^^[{X^) - (X)2]i/2 (7) 



(e.g., Briggs 1993) where X = cos 6, sin^ 6 or cos^ B, and the averages {X) and (A^) are 
taken in the HBM. These error bars, appropriate for comparing observations with the HBM 
(solid lines), are slightly larger than the error bars that would be appropriate for comparing 
with the isotropic model (dashed lines) .|^ 



11 



Briggs et al. (1996) have studied the angular distributions of 1005 GRBs, putting 
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Figure ^ shows that the 3B catalog is statistically consistent with both isotropy and 
the HBM at the 2-o" level or better, for all angular statistics and at all peak photon fluxes. 
Nevertheless, isotropy fits the observations better than the HBM in every graph, unlike in 
our earlier studies based on IB and 2B (LDT; LD96a). At the bright end, the HBM offers 
several clear predictions which may be tested as BATSE accumulates more data. 

4.3. Short &; Long Burst Subclasses in the HBM 

When the numbers of BATSE bursts are plotted as a function of burst duration (i.e., 
a histogram), they show a statistically significant two-peaked distribution (Kouvelioutou 
et al. 1993; Koshut et al. 1996; Meegan et al. 1996). This suggests that GRBs can be 
physically grouped into two classes. In this paper we have analyzed the distribution of 
short and long bursts together, using the photon peak fiux on a 0.256 s time scale, P256, 
as a brightness measure in order to make a homogeneous merged set of data from BATSE 
and PVO. However, in two earlier papers (LDT and LD96a), we studied the BATSE 
duration subclasses independentally, and compared them to HBM predictions. This was 
useful for exploring the possibility that the HBM, or indeed any galactic halo model, might 



constraints on galactic models, including the HBM. They found a deviation of 3.1 cr between 
the HBM and BATSE values of (cos^G) for the brightest 1/8 of their sample, which is larger 
than we find here. Our analysis differs from Briggs et al. in two ways. First, a new intrinsic 
peak luminosity was found from fitting 3B/PV0 (§4.1), which is slightly larger than the 
luminosity that DLT found from fitting IB/PVO, as used by Briggs et al. Second, our 
statistical uncertainties come from equation (0) where the averages (cos^G) and (cos^G) are 
evaluted in the HBM rather than in isotropic models. Isotropic models underestimate the 
relevant uncertainties by more than 10%. 



-25- 

apply to one subclass only. In these studies we used P1024 for the long bursts (Tgo > 2 s) 
and Pq4 for the short bursts (T90 < 2 s). The peak luminosity obtained for the 2B long 
bursts alone was (4.0 ±1.1) x 10^^ photons s^^ (LDT), consistent with the values given in 
eqn. (|]) . Our plots of GRB angular statistics as a function of burst brightness for the two 
subclasses (analogous to Figure 6) showed that the BATSE data was fit by the HBM and 
isotropic (cosniological) models almost equally well. No statistically significant deviations 
from isotropy or from the HBM were found in either subclass at any peak fiux level. 

5. Burster Energy Budgets and Repetition Rates 

If halo GRB sources are observed as soft gamma repeaters (SGRs) when they are 
young, then a rough lower bound on the galactic formation rate of bursters is F > 10~^ 
yr^^. This is based on the fact that 3 SGRs have been detected in likely association with 
supernova remnants of age ~ 1 x 10^ years. Note that SGR activity is episodic (e.g., 
SGR 0526-66 was observed to burst only during the years 1979-1983; Norris et al. 1991), 
and only a small fraction of Type II supernova remnants have confirmed radiopulsar 
(non-SGRQ counterparts (Frail 1997; Frail & Moffett 1993), thus there might exist many 
heretofore undetected SGRs of ages ^ 10"^ yrs in the galaxy. Parameterizing F3 = (F/10~'^ 
yr^-*^), we find: 

0.1 < F3 < 10, (8) 



^^The Prot = 8.0 s rotation period and other properties of SGR 0526-66 clearly indicate 
that SGRs are, as a class, distinct from rotation-powered pulsars. Only one SGR, 1806-20, is 
embedded in a radio plerion, and this has a steeper spectrum than a pulsar plerion (Kulkarni 
et al. 1994). This plerion might be magnetically-powered, as explained in Appendix A. 
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where the upper bound comes from the fact that the rate of burster (SGR) formation 
cannot exceed the total galactic supernova rate. 

If the BATSE sampling depth is D = 100 -Dioo kpc and if all bursters acquire recoils of 
magnitude ~ 1000 V3 km s~^ (cf. Appendix A), then the time for a burster to drift out to 
the BATSE sampling depth is tjj = 10^ Diqq Vf^ yrs. The total number of bursters in the 
volume sampled by BATSE is then 

MD = Ttn = 10' DiooVi'Ts. (9) 

It is worthwhile emphasizing the meaning of this equation. Given the observed SGR 
formation rate [eq. (|^)] and their observed recoils (Appendix A) we can conclude on a purely 
empirical basis that there exists a galactic halo of > 10^ aged SGRs within the distance 100 
kpc (DT92). The number of old SGRs in the halo might be orders of magnitude higher if 
Fa > 0.1, or if a significant fraction of SGRs remain bound to the galaxy (eq. [^ below). 
SGRs are known to emit multiple soft gamma-ray bursts of energy ~ 10^^ ergs when they 
are young; if they continue to emit bursts with similar energies but harder spectra as they 
age, then they constitute halo GRB sources (DT92). 

The (extrapolated) full-sky rate of bursts that are bright enough to be detected by 
BATSE is ~ 800 yr~^ (Meegan et al. 1992). If emitted by aged SGRs, then the mean 
repetition time for bursts from a given star is0 

TanB ^ 110 Ts Vf' D^oo yrs (10) 

This implies that evidence for repetition might be difficult to find in presently available 
data. 



^^ In the HBM, the beaming of bursts reduces the number of observable bursters by only 
a mild factor / ^ 1 - {2/?,)Ro/{<phD) ^ 1 - 0.16 D^^Kcpb/ 20°)-^ ^ 0.84, where Ro ~ 8.5 kpc 
is a galactic disk dimension. This (unimportant) factor is included eq. (10). 
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The intrinsic luminosity of a GRB near the BATSE detection threshold of 
Fbatse — 1 X 10"'^ erg cm"^ s"^ (Meegan et al. 1992) is 

Egrb ^ 1 X 10^1 Dl,, (^) erg, (11) 

where Afi is the solid angle of the burst emissions. This agrees with the luminosity quoted 
in equation (j^), for a best-fit BATSE sampling depth -Dioo = 1-82 and (pb = 20°. The 
effective duration of BATSE GRBs is ~ 10 s, as found from the mean ratio of fiuence to 
peak flux (Lee & Petrosian 1996). Thus the burst energy is 

Egrb ^ 1 x 10^^ ^2^^ (^^^ ^^g. (12) 

In its observable (by BATSE) lifetime of ta = 10^ -Dioo ^~^ yrs, each star emits a 
number of bursts Mb = tol'TGRB = 10^ Fg^^ The total energy required from each burster is 
then Etot = Mb -^grb, or 

E,o, ^ 3 X 10^^ DU V-^ (^-^^) erg. (13) 

This is the basic energy requirement for burst sources in the HBM. Note that in the 
small-angle limit|^ 

/1-COS0A /^y 

V 0.06 / V20°/ ' ^ ^ ' 



^^In the remainder of §5 and in §6.3 we will quote many results with a scaling factor 
(0fe/2O°)^. These equations are given implicitly in the 0^ ^ 1 limit. However, all these 
results can be generalized to any 0^ by using eq. (14). In particular, a magnetar model 
with isotropic (0b = 7r/2) bursters could agree with BATSE and PVO data if the initial 
onset of bursting were delayed (LD92; Bulik & Lamb 1996; Appendix B, eq. [B7]). Energy 
budget constraints for this scenario are given by dropping the (0^/20°)^ factor in eq. (17) and 
multiplying by a numerical factor (0.06)^^ = 16.6. These same isotropic results would apply 
if bursts were beamed, but in random directions; for example, if the alignment scenarios of 
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Rotational energy, the power source for ordinary pulsars, is probably not capable of 
providing this much energy for GRBs (LD92; Hartmann & Narayan 1996). Furthermore, 
there is evidence that SGRs are slow rotators {Prot = 8.0 s for SGR 0526-66). We now 
discuss more promising energy sources. 

5.1. Magnetic Energy 

A neutron star with interior flux density B = 10^^ Biq G, idealized as uniform 
throughout the star, has magnetic energy 

Eb = B^R^/G = 2 X 10^^ ^1^6 erg. (15) 

This is a plausible estimate for the magnetic energy of a young neutron star with an internal 
toroidal field generated by a saturating a-Q dynamo (DT92; TD93). It exceeds by a factor 
~ 10'^ Fa the GRB energy requirement of eq. ([T3|). The transition to superconductivity 
occurs at i? < i?c ~ 3 X 10^^ G. When B < Be, magnetic flux is concentrated by 
supercurrents into discrete flux tubes, and the energy is (PRR): 

Eb = BB^R^/G = 5 x 10^^ 5i5 erg. (16) 

As a magnetar field evolves and decays it will enter the regime of eq. (16), even if it is in 
the eq. (15) regime initially. 

Equations ( |I3D and (|16D together imply that magnetic energy can suffice to produce all 
GRBs in the HBM if the burster internal magnetic field satisfies 

B>6xlO'^DUr^' {^X G. (17) 



§2.2 did not operate. Although the energy per burst would then be reduced by (AQ/Atc), 
one would need {AQ/4tt)~^ more bursts to produce the same rate of observable events; these 
factors cancell in the energy budget (Hartmann & Narayan 1996). 
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This criterion is met, and in fact exceeded by a factor ~ (10-10^) T^, in magnetar formation 
scenarios (DT96) and in the magnetar model of SGRs (TD95, TD96). 

Regardless of the energy source for bursts, if the energy release occurs inside a neutron 
star and is communicated outward by material and Alfven waves (e.g., Blaes et al. 1989), 
then very strong magnetic fields, B > Bq are probably necessary. Weaker fields would 
allow vibrational crumbling of a halo neutron star's crust and thereby produce an X-ray 
excess (Appendix B). 

5.2. Accretion Energy 

The sudden accretion of a planetesimal (asteroid or comet) with mass 

M, = 5 X 10^1 Dl, e-' (^) gm (18) 

onto a neutron star would release the energy of a halo GRB, where e^ is the efficiency of 
gamma ray production. This is a scaled-up version of the GRB theory of Harwit & Salpeter 
(1973) and Colgate & Petchek (1981). 

To power all 10® T^^ bursts in a burster's observable lifetime, the total mass of accreted 
planetesimals is 

Mtot = 5 X 10^^ Tg-i Dl, e-' (^) gm, (19) 

similar to the mass of the Earth. These planetesimals presumably condense from a 
gaseous disk. Such a disk (M^isk > ^tot) could be captured by a HVNS when it overtakes 
clumpy supernova ejecta (Woosley & Herant 1996) or experiences a recoil directed toward 
a stellar companion (Colgate & Leonard 1994). Alternatively, a disk could form when 
a rapidly- rotating proto-neutron star sheds excess angular momentum (e.g., Thompson 
& Duncan 1994). In this case, the intial radius of the retained disk could not exceed 
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-Rj,disk ~ {GM^jV^) ~ 2 X 10^ V^ ^ km; but it would subsequently experience viscous 
spreading. 

It is not obvious how the beaming and alignment requirements of the HBM are met 
in such an accretion-powered scenario, although the physics is complex enough that a 
mechanism might be identified. Alternatively, GRBs could be unbeamed, or beamed in 
random directions, but with a delayed onset of bursting (LD92). Colgate & Leonard (1994) 
estimated that a time interval ~ 10^ years, as required to fit BATSE data, occurs before 
planetesimals condense in the disk. In the pulsar-planet formation scenario of Lin, Woosley 
& Bodenheimer (1991), some objects as large as ~ 10^^ gm formed much earlier, after ~ 10^ 
yrs, but ^10^ yrs were needed to form planets. Planets might be necessary to defiect a 
significant number of objects toward the neutron star (e.g., van Buren 1981; Tremaine & 
Zytkow 1986), so a delayed turn-on as Colgate & Leonard propose seems feasible. 

Although the gross energy requirements of halo GRBs can easily be met in an 
accretion-powered model, these models have two potential problems: (1) producing hard 
non-thermal gamma ray spectra without over-producing thermal X-rays; and (2) producing 
a large enough rate of abrupt planetesimal captures (eq. p^). 

Some form of magnetically-mediated gamma-ray production is usually invoked to 
solve the first problem (e.g., Colgate 1992). That is, the accreting matter distorts the 
magnetic field, which — via nonlinear Alfven wave damping, tearing instabilities, or some 
other mechanism — accelerates charged particles to produce gamma rays. Note, however, 
that the total exterior magnetic energy of a neutron star with -Bdipoie = 10^^ B12 G is less 
than the energy of a single GRB in the galactic halo by 

^ = 02BI,D-SA^)\ (20) 

J^GRB V 47r / 

It is difficult to understand how such a relatively large accretion energy (eq. [jT^) could 
be channeled through such a weak field without enormous inefficiency, resulting in the 
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over-production of thermal X-rays. This problem is discussed in more detail in §7.3.1 of 
TD95.|3 

A strong magnetic field is also generally invoked to increase the capture cross-section 
of a neutron star, thereby enhancing the putative GRB rate (Colgate & Petchek 1981; van 
Buren 1981). A crude estimate of the radius Re at which magnetic capture will occur in a 
single orbital pass is given by the condition that the kinetic energy of mass Mp at Re be 
comparable to the magnetic dipole field energy external to Re- This yields 



4 d2 X 1/2 / , X -1 

1.7 B,,D-,^,e\'^tA (21) 



Re _ ( R* -^dipole \ _ 1 7 D n-1 .1/2 ( '^b 

r: - [ggmJTJ - ^-^ ^'' ^^°° '^ 20 



'*^"p, 



for i?^ = 10 km and M^ = 1.4 Mq. Thus the single-orbit capture cross-section is barely 
enhanced over the geometrical cross-section of the star for pulsar-strength fields and halo 
locations. Multi-orbit magnetic captures also occur (van Buren 1981); but as shown by 
Woosley & Herant (1996), very oblique captures produce extended soft X-ray events. 

For local galactic disk bursters (-Dioo ~ 10^'^), the context in which planetesimal 
accretion was originally proposed, neither of these problems (as roughly illustrated by the 
small numbers in eq. |^ and |^) are serious. However, for halo GRBs, dramatic beaming 



{AQ <^ 1) or magnetic fields much stronger than 10^^ G, or both, seem theoretically 
desirable. 

Thus we suggest that planetesimal accretion onto magnetars may power halo GRBs. 
Some subset of the rapidly-rotating proto-neutron stars which become magnetars (according 
to the conjecture described in DT96) may shed excess angular momentum into low-mass 
disks (eq. [l^), or capture disk material as they move through clumpy supernova ejecta 



(Woosley & Herant 1996). The exterior magnetic energy of these stars is greater than the 



^^TD95 focus on accretion-powered models of the 1979 March 5th event's hard initial 
spike, but the spectral requirements are similar (FKL). 
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energy of a halo GRB by ~ 2 x 10^ Bf^ D^qq (Af2/47r)^^, potentially catalyzing substantial 
non-thermal emissions. The single-orbit capture cross-section for planetesimals is enhanced 
over the geometrical cross section by a similar factor. 

The HBM could apply in a magnetar-accretion model because the mechanisms for 
stellar recoils and for aligned beaming described in §2 would hold. Alternatively, a time 
delay for planetesimal formation (Colgate & Leonard 1994) may operate (LD92). This 
physical model for halo GRBs is viable even if all magnetically-powered activity terminates 
after the SGR phase, as long as a magnetar retains a substantial fraction of its initial dipole 
field, anchored in the stably-stratified core, for a time ^ 10® V^^ -Dioo years. 

6. Model Extensions and Refinements 

6.1. Intensity Variations Across the Beam 

So far we have assumed that the gamma ray intensity is uniform within a beaming cone 
of half opening-angle (pi,. In this section we will study the implications of a gamma-ray 
intensity which varies as a function of the angle from the beam axis. 

We will use the relativistic beam model of Mao & Yi (1994). In this model, gamma-ray 
emitting material is expanding with Lorentz factor 7 = (1 — t>^/c^)~^/^ into a cone with half 
opening angle A^. Emissions are assumed to be isotropic in the rest-frame of the emitting 
material, which means that they are beamed to within an opening angle ~ 7"^ about the 
material velocity vector in the observer's frame. Integrating over the material beam and 
over photon emission angles, Mao & Yi (1994) derived the net observable photon intensity 
as a function of angle 6 away from the beam axis: 



27r rA8 

d(t)' 

Jo 



m = A / rf0' / sin^' dQ' X [7(1 - [vx/c])]-'^^^'^ . (22) 
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In this equation, x — cos^' cos0 — sin6'' sin0' sin^; 5 is the photon spectral index (we will 
take (5 = 2); and A is a normalization factor, which we will define such that /(O) = 1 for all 
values of A6' and 7. 

This idealized physical model provides a useful two-parmeter family of beams with 
centrally-peaked brightness. Yi and Mao (1994) used this model to show that relativistic 
beaming effects in galactic halo GRB models can produce anti-correlations between the 
brightness and duration of observed bursts, mimicking cosmological time dilation.]^ Such 
correlations were found by Norris et al. (1993) in the BATSE data, but not by Mitrofanov 
et al. (1996). 

In Figure ^ we plot I {9), as given by equation (^2]). We have kept A^ = 20° in this 
plot, which ensures that the beam half-width at half-maximum stays close to 20°. (This 
is likely to fit BATSE and PVO data best; see eq. [3].) Figure 7 shows that, for small 7, 
bursters become potentially detectable at large 6 (> 20°), although only at low intensity, 
i.e., only within a diminished sampling depth by BATSE. The 0{, = 20° uniform beam 
case, considered in §3 and §4 (solid line in Figure 7), is equivalent to the limit 7 ^ 00 on 
this plot. 

What values of 7 are plausible for a halo GRB model? A pure photon/pair fireball 
(presumably channeled by the magnetic field into a restricted solid angle) with intrinsic 



luminosity appropriate for a galactic halo GRB would have 7 ~ 10 or slightly less (|Piran 



|fc Shemi 1993 ). This is comparable to the value used by Yi and Mao (1994) in their halo 



"time dilation" study. The Lorentz factor could be moderately reduced, say to 7 ~ 3, by a 



^^When bursters are unbound from the Galaxy, as in the HBM, more distant sources are 
systematically older. Any trend of spectral softening or increasing burst duration with age 
could then mimic cosmological time dilation. 



-34- 

plausible contamination of baryons.[] 

To study the implications of these I{6) models, we assume as before that the axi 
of beamed GRB emissions are aligned with peculiar recoil velocities Vr in a galactic 
population of bursting stars, and we do new Monte Carlo calculations of the observable 
GRB distributions. Specifically, for all the I {6) functions shown in Figure 7, we repeated 
the log A/" — log P model fitting against PVO and 3B data (as shown in Figure 4 for 
the uniform beam case), finding the best-fit peak photon luminosity, now defined at the 
beam center (6* = 0). We then used this best-fit value to calculate the cumulative angular 
statistics as a function of peak photon flux, similar to Figure ^. We show these results in 
Figure ^. For clarity, we have plotted the difference of the angular statistics (cosG), (sin^6) 
and (cos^B) from the values predicted in the uniform beam case (i.e., the (pb = 20° HBM 
values, as plotted in Figure |^), for various intinsic beaming patterns I{0). 

Figure 8 shows that, as long as the half- width at half- maximum of the gamma-ray 
beam is held constant, observable results for various beam intensity patterns I{6) deviate 
little from the uniform beam case. (Note that the vertical axis scale in Figure ^is expanded 
by a factor ~ 10 relative to that of Figure |^.) This is our main conclusion. In fact, all 
of the I{6) beaming patterns considered in Figure 7 and Figure 8 fit the present 3B data 
with better than 3-cr confidence at all observed peak photon flux levels, as can be seen 
by comparing Figure 8 with the BATSE data of Figure |^. Thus a uniform beam is not 
required in the HBM, as long as the variation within the beam is not too large. 



^^ Of course, fireballs produce (blueshifted) thermal-spectrum photons, which is not what 
GRBs look like; however this gives at least a first estimate of what the hydrodynamically- 
driven Lorentz factor would be. Electrodynamic acceleration (e.g., driven by magnetic 
reconnection or Alfven waves) might produce larger 7, and/or different beam patterns. 
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6.2. Burst Luminosity Functions 

We now relax the "standard candle" assumption (point [4] in the second paragraph of 
§2.1). We will consider GRB luminosity functions of the form 

$(L) OC L-^ for Linin <L< L^nax- (23) 

Constraints on GRB luminosity functions have been previously inferred from BATSE and 
PVO data by Ulmer, Wijers & Fenimore (1995, hereafter UWF), Hakkila et al. (1995), and 
by other workers cited therein. In these studies, the bursters' spatial density distribution in 
the galactic halo was assumed to be 

The intensity distribution of GRBs (log M - log P curve) in the HBM resembles that of 
the a = 2 case of this model, although the angular distribution of observable bursters as a 
function of r is very different. 

These previous studies concluded that the acceptable dynamic range (i^max/-^min) for 
GRBs observed by BATSE is quite limited. In particular, various studies agree that 90% of 
detected bursts must have peak luminosities within a factor of ~ 10, for a wide range of 
indices a and 13. The basic reason for this is that the turnover of the observed log M - 
log P curve (Fig. 4) from a slope d\o^J\f {> L) / d\ogL = —3/2 (the "Euclidean core") to a 
shallower slope at fainter fluxes happens rather sharply. A GRB luminosity function that 
spanned much more than one order of magnitude would smooth out this turnover to an 
extent inconsistent with observations. 

In the a = 2 case (relevant to the HBM) UWF found, furthermore, that the dynamic 
range of the luminosity function, i^max/-^mini can only be as large as ~ 10, and still flt 
observations, if the power-law index lies in the range (3 = 2 — 3. (See Fig. 2 in UWF). 
Other values of (3 would require signiflcantly more restricted dynamic ranges. This result 
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can be understood analytically. If we assume L^;^^/ Lq = Lo/Lmin = i^ where Lq is the best 
standard-candle (mono-luminosity) fit to the data, then simple integrals show that the 
average luminosity satisfies L/Lq ^ 1 for 2 < k < 20 only when /3 ~ 2 — 3. L/Lq is very 
different from unity when k > 3 for all other values of j3. Since luminosity functions with 
L/Lq ~ 1 give the best fits to observations, the dynamic range, given by k"^ = LmaxZ-^min, 
can only be large for 2 < /? < 3. More general analytic constraints on GRB luminosity 
functions have been derived by Horak et al. (1996). 

These above-quoted results were all obtained from fits to the the observed log Af — log P 
curve only. Next we consider constraints on the GRB luminosity function from the observed 
log A/" — log P distribution combined with the observed angular distributions. These 
combined constraints were previously considered by Hakkila et al. (1995) for a galactic halo 
of bursters with density given by eqn. ( ^4]) and no preferred beaming directions. Here we 
focus on the very different case of the HEM. (Of course, the GRB angular distribution 
cannot constrain the luminosity function of cosmological GRBs, since the bursters are 
distributed essentially isotropically around Earth at all intensity levels.) 

Intuitively, luminosity functions could either wash out or enhance the angular 
anisotropics predicted by HBM for bright subsets of GRBs (Fig. 6). If the luminosity 
function of bursts is a positive (negative) power-law, then at a given observed flux level, a 
greater (smaller) number of more distant bursters are counted than in the mono-luminosity 
case. Consequently, the degree of anisotropy is reduced (enhanced) since the HVNS 
distribution appears more and more isotropic at greater distances from the Earth. 

To illustrate this, we consider two cases: (/5, k) = (—3,2), (3,4.4). Both models 
give acceptable fits to the log Af — log P distribution of the combined PVO and BATSE 
data, but the allowed dynamic range L^ax/Lmm = ^^ is much smaller in the /3 = — 3 case, 
for reasons discussed above. Figure shows corresponding angular distribution statistics 
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(dashed and dotted curves) compared to the mono-luminosity case (thick sohd curve). It is 
clear that the luminosity function with positive slope {—f3 = 3) affords a better fit to the 
angular position data than the standard-candle model does. For j3 = 3, on the other hand, 
the angular statistics deviate from BATSE by ~ 3-(T at the faint end. (Plotted error bars 
are 3-cr.) 

We conclude that only a limited spread of intinsic luminosities among GRBs observed 
by BATSE and PVO are admissible in the HBM, as in other galactic halo models. 
Depending on the exact form of the luminosity function, this could either improve or worsen 
the model's fit to the BATSE angular distribution. However, in either case the effect is 
limited because the dynamic range is small. Equivalently, one could say that in a fit to 
the combined GRB brightness and angular position data, angular information increases 
(decreases) the statistically-allowed dynamic range of the luminosity function when (3 < 
(/9 > 0), but this dynamic range is less than ~ 10 in any case. 

Note that we have considered only one possible generalization of the basic (standard) 
HBM here: a power-law luminosity function that is uniform among bursters. Other 
conceivable generalizations, such as luminosity functions that vary with burster age, or 
luminosities correlated with 0{, or Vr, would produce different results. 

6.3. Weakly-bound Bursters in the HBM 

Next we consider relaxing the assumption that all bursters are born with Vr = 1000 
km s~^ (point [2] in the second paragraph of §2.1). If some fraction of bursters are born 
with velocities less than the local galactic escape velocity, which is about Vesc ~ 600 km s~^ 
at typical birthsites in the disk, these stars will orbit the galaxy. Podsiadlowski, Rees and 
Ruderman (1995) have emphasized the important point that these weakly-bound bursters 
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do not return to the small galactic radii at which they were born because the gravitational 
potential in the halo is probably markedly non-spherical. Large deviations from spherical 
symmetry are predicted by almost all models of galactic halo formation (e.g., Dubinski 1992 
and references therein). Thus weakly-bound bursters (i.e., bursters born with velocities 
narrowly below the local galactic escape velocity) will tend to stay in the outer halo, 
following wide-ranging, non-periodic orbits. After one or more orbits, their magnetic axi 
orientations are randomized. If these weakly-bound bursters continue to emit (beamed) 
bursts, their contribution to observed GRBs must be added to the bursts from the escaping 
ones (i.e., the unbound "standard HBM" sources, studied in §3 and §4). This could, under 
some conditions that we will discuss below, produce a net angular distribution of GRBs 
that is more isotropic than that of the standard HBM, especially at the faint end. However 
the basic HBM mechanism for producing the "Euclidean core" in the ORB brightness 
distribution, as observed by PVO, may still apply. 

To quantify this discussion, assume that a fraction / = 10^^/_i of all bursters are born 
with velocities below the local galactic escape velocity, and the remainder have Vr = 1000 V3 
km s~^ as before. This is a crude way to estimate the effects of a realistic (Gaussian) 
velocity distribution. If the total birthrate of GRB sources in the galaxy is F = 10^^ Fa 
yr^^ (eqn. [Q), then the number of bound bursters presently orbiting in the halo is 

ATbound = F / Tgal = 1 X 10^ F3 /-I, (25) 

where Tgai ~ 10^° yr is the age of the galaxy (i.e., the period during which bursters have 
been forming). The relative abundance of bound bursters compared to unbound ones 
(eqn. |Q) within the BATSE sampling depth oi D = lOODioo kpc is then 

10 /_i Dfoo ^3. (26) 



■'Abound 



■'V unbound 

Thus if > 1% of bursters remain bound to the galaxy (/ > 10^^), bound stars will 
numerically dominate the halo (DT92). If the bursters continue to emit bursts within 
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beaming cones of half angle (pb, but with orientations randomized by their orbits, for a 
bursting lifetime rufe = 10^° tiq years, then the fraction of all bound bursters that are 
observable at Earth is 

Multiplying equation ( pBf ) and equation (^) yields the relative number of observable bound 
bursters compared to unbound ones: 

^6 = 0.6 /_! no Droo ^3 (^) . (28) 

It is interesting that this can be of order unity. 

A fraction ^^^/(l + "^b) of all observed GRBs will come from bound (orbiting) bursters. 
How anisotropic will be the sky distribution of these "bound bursts" ? This depends upon 
several uncertain factors, especially the burster recoil velocity distribution and the form of 
the galactic halo potential. Here we limit ourselves to a few general comments, based upon 
the results of PRR. 

If many bursters are born with Vr ^ 500 km s~^, significantly below the escape velocity 
from the center of the galaxy, then the bound burst distribution will be very anisotropic 
(cf. Fig. 2a and 2b of PRR). If, however, most bound bursters aquire velocities only 
narrowly below the local escape velocity, then the bound burster distribution could be 
nearly isotropic, as shown in Fig. 5b of PRR. This figure shows the burster distribution for 
kick velocities uniformly distributed on the range 700 < Vr < 900 km s~^, and a realistic, 
nonspherical galactic potential. 

In the magnetar model of GRBs, burster recoils are imparted by a complicated series of 
magnetically-induced neutrino impulses (§2.2), with durations comparable to the coherence 
time of magnetic structures near and below the neutrinosphere (DT92; TD93). Because 
these impulses are additive, the central limit theorem implies that a Gaussian with mean 
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V and variance ay is probably a good first approximation to the Vr distribution among 
bursters. The precise values of V and av are uncertain, although SGR data suggest 

V ^ 1000 km s^^ (Appendix A). Thus the Vr distribution for bound bursters is probably 
the low-Vr tail of a Gaussian. If ay ~ 200-300 km s~^, then the velocity distribution of 
bound bursters is sharply peaked just below the escape velocity, and conditions for a nearly 
isotropic sky distribution of bound bursts could be met. The distribution of bound bursters 
might still show mild "hot spots" in some presently-unknown directions associated with 
distortions in the outer galactic halo potential (PRR), but these anisotropics might not be 
apparent from the statistical tests used in this paper. [More general statistics (e.g., Briggs 
et al. 1996; Tegmark et al. 1996) would be useful for constraining this effect.] 

Assuming that ay ~ 200-300 km s~^ and that other conditions for a nearly-isotropic 
bound burst distribution are met, then the total model galactocentric quadrupole moment 
would be 

(cos'e) - 1/3 = (^b + l)-i [(cos2e)HBM - (1/3)(1 - ^b)], (29) 

where (cos^9)hbm is the moment in the "standard HBM", as shown in Figure (6), and we 
have assumed that the bound burster distribution has (cos^0) = 1/3. Thus for \E'(, = 1, 
the deviation from isotropy would be reduced by (1/2) compared to the standard HBM. 
Similar results might apply for (cosB) and (sin^6); however it is also possible that cTy ^ 300 
km s~^ and that the addition of bound stars makes the model fits to BATSE observations 
worse. We will study this issue in a future paper. 

Note that the mean time between bursts from a given star is longer by a factor (1 + "^f,) 
when a population of bound bursters is present: 

TauB = 110 (1 + ^b) Ts ^3^' I^ioo yrs (30) 



cf. eq. ll^])- Despite this fact, the energy required for each burster is significantly greater, 
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because of the long bursting lifetime: 

Et,t = 5 X 10^« rio I^ioo Tg"^ V^ (1 + vl/,)-i (^\ erg. (31) 



This still might be supplied by magnetic energy (cf. eq. [jT5[). Alternatively, if the interior 
flux density satisfies -B < 3 x 10^^ G initially, then from eq. (|16D and eq. (|3T|) we infer that 
magnetically-powered bursting would have to terminate after an active lifetime no longer 
than 

T^,, = W' 5i5 Ts 1^3^^ ^roo (1 + ^b) [0) yrs. (32) 

This corresponds to tiq = 0.1 in eq. (|^), suggesting that magnetically-powered, orbiting 
bursters might cease GRB activity too early to make a substantial contribution to the 
observable burst distribution {'^b ^ !)• However, accretion-powered bursters (§5.2) would 
not be subject to this bursting lifetime limit. 



7. Concluding Discussion 

7.1. Summary 

In this paper we have considered models for GRBs from high-velocity neutron stars in 
the galactic halo. In particular, we compared GRB observations with a model having the 
following properties: [1] bursters are born at positions distributed like young Pop. I stars 
in the galactic disk; [2] with randomly directed recoils Vr = 1000 V3 km s~^; [3] they emit 
GRBs at a constant rate, with [4] constant luminosity, and [5] the gamma ray emission is 
beamed parallel and anti-parallel to Vr, within an angular radius 0^. 

This "standard halo beaming model" (S-HBM) adequately fits all published data on the 
distributions of GRBs from the BATSE and PVO experiments, when 12° V^^ < (pb < 30°. 
Models for GRBs from sources at cosmological distances fit present BATSE observations 
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better than the S-HBM does, but not by a statistically-compelhng margin. These fits to the 
data, for Vr = 1000 km s~^ and (ph = 20°, are shown in Figure 4 and Figure 6. Note that 
the fit improves for Vr > 1000 km s~^ (as found in most SGRs; Appendix A). 

An important feature of the HBM is that almost all bursters born in the Andromeda 
Galaxy (M31) or other nearby external galaxies are unobservable at Earth because of 
misdirected beaming (LDT; Li, Fenimore and Liang 1996). Only when the samphng depth 
of the instrument reaches M31 can the effects of M31 possibly become evident. This means 
that the often-quoted upper bounds on the BATSE sampling depth based on considerations 
of M31 (e.g., LD92; Hakkila et al. 1994) do not apply to the HBM. Bounds on halo models 
based on the nondetection of X-ray bursts from M31 (Li & Liang 1992) or in samples of 
nearby external galaxies also do not constrain the HBM, as long as the X-rays associated 
with GRBs are beamed in a similar way as the gamma rays. 

We have also considered three refinements to the S-HBM, relaxing the assumptions 
[2], [4] and [5] in the above list: 

Intensity variations across the beam. We analyzed models with centrally-peaked 
gamma-ray emission beams; i.e., with intensity declining as a function of angular separation 
from the beam axis, as in relativistic emission models. We found that the observable model 
results are not very sensitive to the details of the beam intensity distribution. Adequate 
fits to the data were obtained as long as the angular radius at which the intensity falls to 
half its central value is ~ 20°, and the angle at which it drops by another factor ~ 2 is not 
many times larger. 

Burst luminosity functions. We explored the consequences of a spread in GRB intrinsic 
brightnesses, considering luminosity functions of the form $(L) oc L^^ for Lmin < L < Lmax- 
We found that, in the HBM as in other galactic halo models, only a limited dynamic range 
-^max/-^min < 10 or 20 is admissablc among GRBs observed by BATSE and PVO. Depending 
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on whether the power-law index satisfies (3 < or (3 > 0, the spread in luminosities either 
improves or worsens the model's fit to the BATSE angular distribution. However, in either 
case the effect is small because the dynamic range is limited. 

Models with a subpopulation of bound bursters. In the S-HBM, all bursters are born 
with recoil velocities Vr = 1000 V3 km s~^. Realistically, one expects bursters to aquire a 
distribution of K values. The low-K tail in this distribution could correspond to bursters 
that remain weakly-bound to the galaxy. Podsiadlowski, Rees & Ruderman (1995) have 
shown that these bound stars follow non-periodic orbits in the non-spherical galactic halo 
potential, generally not returning to the small galactic radii where they were born. After 
one or more orbits, their magnetic pole (beaming axis) orientations are randomized. If they 
continue to emit (beamed) GRBs at late ages, then the bound and escaping bursters might 
produce comparable numbers of observable GRBs. (The exact ratio, eq. p8|], depends upon 



several uncertain model parameters.) The observable distribution of GRBs would then be 
an appropriately weighted sum of the S-HBM burst distribution, describing the unbound 
bursters (Fig. 4 and Fig. 6), and the bound burster distribution studied by PRR. The 
resultant model angular distribution might be significantly more isotropic than that of the 
standard HBM, if the orbiting bursters are only weakly bound to the galaxy (see §6.2). 
This plausible refinement of the standard HBM merits further study. 

7.2. Conclusion: Cosmological or Galactic GRBs? 

Our main conclusion is as follows. The models studied in this paper give acceptable 
fits to observations as long as certain special conditions (e.g., uniform bursting rate, 
12° Vf^ < (pb < 30°) are satisfied. These conditions, although possible and arguably 
plausible, have no clear and compelling theoretical justification. Other special conditions of 
varying plausibility are required in other halo models (e.g., LD92; see Hartmann et al. 1994 
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for a comprehensive review of halo models). The cosmological hypothesis, on the other 
hand, fits present observations generically, and thus is much more attractive in this sense. 
Moreover present data are generally better fitted by the cosmological models than by the 
S-HBM, as shown by Fig. 6. 

Based upon this powerful argument, we believe cosmological GRB models to be more 
promising than galactic halo models, given the evidence available at the time of writing. 
However, some caution is appropriate because there exists no compelling theoretical 
understanding of GRBs at either distance scale. The study of GRBs has taken many 
unexpected twists and turns in the past: only five years ago there was a wide consensus 
favoring neutron stars at distances ^ 10^ pc. 

If the cosmological hypothesis is correct, then continued BATSE data collection should 
be able to rule out the standard HBM using the test of Fig. 6. A doubling of the data base 
beyond the 3B catalog might rule out the S-HBM with > 3-cr confidence. 

A plausible refinement of the standard HBM that might be more difficult to falsify 
includes some bursters that are weakly-bound to the galaxy (§6.3). Such a model can be 
jointly constrained by limits on GRBs — or associated X-ray bursts — from nearby external 
galaxies (e.g., Li & Liang 1992), which limit the bound burster population, and BATSE 
angular position observations, which more strongly limit the unbound ones. We will 
quantify these constraints in a future study. 

If GRB do come from the galactic halo, then deviations from isotropy on the sky will 
eventually manifest themselves. In any case, continued operation of BATSE is essential for 
making these tests. Checking BATSE with an independent experiment is also desirable. It 
is important to verify that systematic errors do not wash out the very low-level anisotropies 
predicted by halo models. 
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GRB repetitions, if verified, could lend credence to halo models; however the predicted 
mean time interval between repeat bursts from old SGRs in the halo is Tq^b ~ 100 years, 
with some dependence on uncertain model parameters (eq. [|l^ and eq. 0]), so repeat 
events may be difficult to detect unless they are strongly clustered in time.|3 [^ 

The durations of BATSE bursts were found to be anti-correlated with their intensities 
by Norris et al. (1993) [but see Mitrofanov et al. 1996]. This could be due to cosmological 
time dilation. However, this correlation, if verified, does not give strong evidence against the 
HEM. In the HBM more distant burst sources are systematically older, so any significant 
trend of increasing burst duration or spectral softening with age could mimic cosmological 
time dilation. 

What are the strongest arguments in favor of galactic halo GRB sources? There is 
strong empirical evidence that a large number of old soft gamma repeaters exist in the 
extended galactic halo (DT92). We estimate the number of aged SGRs within 100 kpc to 
be between ~ 10^ and ~ 10^, depending on values of the scaling parameters in equations 



^'^Magnetically-powered GRBs would tend to cluster in time, like stellar fiares, so incidents 
of GRB repetition should eventually be observed if the bursters are isolated magnetars. This 
could be an important observational diagnostic, because most cosmological GRB models do 
not allow burst repetitions. 

^^ Many bursts in the BATSE catalog show distinct episodes of bursting activity separated 
by ~ 10-100 s intervals containing no detected emissions (Lingenfelter, Wang & Higdon 
1994) . It would be interesting to search for such recurrences on somewhat longer timescales, 
~ 10^ s to 10^ s, in order to place bounds on how long the time window for recurrence lasts. 
Unfortunately, BATSE is biased against recording data during a readout period lasting from 
~ 100 s to ~ 90 minutes following any burst, when it it downloading data. A large fraction of 
bursters could produce lower-intensity recurrences during this period without any detections. 
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(D and (pGl). These halo stars have many properties required of classic GRB sources, 
as first pointed out by DT92. They are clearly distinct from radio pulsars (footnote 12 
and Appendix A), thus avoiding the problem of over-concentration of bursters toward the 
galactic disk which plagues many theories based on neutron stars. Moreover, of all known 
astrophysical objects, SGRs have produced the emissions which most closely resemble 
GRBs. They have produced at least one burst of photons with a classic GRB-like spectrum 
(on March 5, 1979; FKL); and when they are young (ages ^ 10^ yrs), they emit large 
numbers of softer-spectrum bursts with energies ~ 10^^ ergs, comparable to the energies of 
GRBs in the halo. 

We have briefly outlined some candidate physical models for these bursters (§2.3, §2.4, 
and §5). One promising class of models involves neutron stars with ultra-strong magnetic 
fields B > 10^^ G, much stronger than the fields known to exist in ordinary radiopulsar 
magnetospheres. Such strong magnetic fields could provide the energy for fiare-like halo 
GRBs (§5.1), or they could mediate the conversion of accretion energy into gamma rays 
(§5.2). Fields this strong could also drive crust fracture events in young neutron stars, 
producing SGR activity (TD95). [Note that the distributions of SGR burst energies and 
waiting times, and their time-correlations, are consistent with crust fractures; Cheng 
et al. 1996.] Finally, strong magnetic fields can drive neutrino emission anisotropics in 
newborn neutron stars, imparting recoils that propell the stars into the galactic halo (§2.4). 
Note that all of these mechanisms could conceivably operate in a non-HBM context, i.e., a 
halo GRB model in which the bursts are unbeamed or beamed in random directions (see 
footnote 14). 

If aged SGRs do not emit classic GRBs, then it is still likely that these stars populate 
the galactic halo (DT92). If they are isolated magnetars, then they will rotate very slowly 
after leaving the galactic disk. They will cool rapidly once their magnetic activity ceases. 
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becoming dark halo relics that are difficult to detect. 

Several promising models of cosmological GRBs also involve ultra-strong magnetic 
fields (Thompson 1994 and references therein). As a general result, gravity-powered MHD 
flows in compact neutron objects can create fields of strength B'^/Sn ~ eGM'^/R^, where e 
is the efficiency. This implies 

B.lO-f^)"'^ m (^y G. (33) 

Although strong by most standards, this is vastly below the fundamental limit of ~ 10^^-10^^ 
G at which magnetic fields become unstable to the vacuum creation of GUT-scale or 
Planck-scale monopoles. 

Thus, the ~ 10^^ G dipole fields of radiopulsars are a million times weaker than the 
strongest possible flux densities attainable, however briefly, in neutron star formation and 
coelescence. The mystery of gamma-ray bursts may require us to flnally recognize and 
understand this new regime of strong magnetism. 
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A. Observations of Soft Gamma Repeaters and Evidence for Recoils 

We have suggested that old SGRs in our galaxy are the sources of classic GRBs. 
This requires that SGRs aquire sufficiently large recoils to propell them into the galactic 
halo. Here we briefly review the evidence for large recoils in SGRs, and discuss the related 
issue of whether SGRs have normal stellar companions. We also briefly discuss how SGR 
observations might be accounted for in the magnetar model (which predicts large recoils, 
§2.4). 

SGR 0526-66 The error box for the 1979 March 5 event QCline et al. 1982|) and the 
SGR-like bursts which followed in its aftermath (Golenetskii, Ilyinski & Mazets 1984) lie 
within supernova remnant N49 in the Large Magellanic Cloud. A point-like X-ray sourceQ 
inside this error box ( Rothschild, Kulkarni fc Lingenfelter 199"^ ) may have a thermal 



spectrum, since no hard spectral component has been detected (Murakami 1996; Marsden 
et al. 1996). Its intensity is consistent with thermal emission from a magnetically-heated 
neutron star (TD96). If the SGR was born at the center of the supernova that produced 
N49, then the displacement of the X-ray source from the remnant center implies a transverse 
velocity (DT92) 

V± = 1200 ± 240 km s-\ (Al) 

where we have used the estimated age of N49, tN49 ~ 5400 yrs (Vancura et al. 1992). The 
statistically expected 3-D recoil velocity is larger by a factor (3/2)^/^, or l^ = 1500 km s~^. 

No radio, infrared or optical emissions have been detected from SGR0526-66 down to 
faint limits (Dickel et al. 1995; Marsden et al. 1996). This is not surprizing since a neutron 
star born with the velocity of eq. (Al) could not remain bound in a binary (DT92); and 



^° That is, a source with an observed X-ray image that is statistically consistent with the 
instrumental point-spread function. 



-49- 



a neutron star with an 8.0 second rotation period, as clearly shown by the 1979 March 5 
event (Mazets et al. 1979), could not produce a significant rotation-powered plerion. 

SGR 1900+14 This burster has been localized to a position lying just outside the 
edge of the galactic supernova remnant G42.8+0.6 (Hurley et al. 1994). Like SGR 0526-66, 
it has a candidate point-like X-ray counterpart, just outside the burst error box (Hurley 
et al. 1996). Assuming that the SGR was at the center of the supernova that produced 
remnant G42.8+0.6, one finds (Hurley et al. 1996) 

where D is the distance from Earth, and tg^p is the age of the supernova remnant (SNR). 
This implies V ~ 2100 km s^^. 

A double infrared source, probably a bound pair of M5 supergiant stars which is 
heavily reddened by interstellar dust, has been suggested as a counterpart object to SGR 
1900+14 (Vrba et al. 1996). This object lies 18" from the center of the observed X-ray 
image (Hurley et al. 1996) and 50" from the burst error box (Vrba et al. 1996) in a region 
of dense galactic starfields, so the identification is not certain. If the infrared binary is a 
true counterpart to the SGR, then the association of the SGR with the supernova remnant 
is almost certainly spurious, because no known process could impart a recoil of magnitude 
(A2) to a pair of M5 supergiant stars. However, the other SGRs are unambiguously 
associated with ~ 10^-year-old SNRs. Thus we favor the interpretation that the infrared 
binary is not physically associated, and the SNR is. In this case, SGR 1900+14 closely 
resembles SGR 0526-66 in all its known (observed) properties. 

SGR 1806-20 This burster is located within the ~ 10^-year-old SNR GlO.0-0.3 
(Kulkarini et al. 1994). It also has an X-ray counterpart (Murakami et al. 1994) which is 
point-like, but with a measured non-thermal spectral component (Sonobe et al. 1994), unlike 
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the other SGRs. Furthermore, there is a synchrotron radio plerion surrounding the SGR 
(Kulkarni et al. 1994) indicating that the burster drives an outflow of relativistic particles. 
The radio spectral index of the plerion is steeper than that of known rotation-powered 
pulsar plerions (Kulkarni et al. 1994). 

In the magnetar model, this plerion could be powered by SGR bursts, as suggested 
by Kulkarni et al. (1994) and studied by Tavani (1994) and Harding (1996). Another 
possibility is that the radio nebula and non-thermal X-rays are powered by quasi-continuous 
low-amplitude Alfven wave emissions from magnetically-induced deep crustal fractures 
(TD95). We favor this second interpretation, for reasons explained in TD96. The other 
SGRs do not power radio plerions; this might mean that the deep crustal field in other 
SGRs is strong enough to overwhelm lattice stresses, driving plastic creep rather than 
quasi-steady low-amplitude seismic activity (TD96). 

Thus in our favored interpretation, SGR 1806-20 has the weakest magnetic flux density 
in the deep crust of known SGRs. Alternatively, it could have the most spatially disordered 
magnetic field in the deep crust, which would also drive a higher rate of crustal fractures. 
In either case, we would expect it to also have the lowest recoil velocity according to the 
mechanisms of §2.4. The displacement of SGR 1806-20 from the center of the SNR is 
difficult to estimate because of the irregular shape of the remnant and because of late 
energy injection into the plerion, but Kulkarni et al. (1994) roughly estimate V± ^ 500 km 

This velocity is small enough, and uncertain enough, that the burster could have 
remained bound to a companion if it formed in a tight and massive binary. In fact, a massive 
companion star has been detected (Kulkarni et al. 1995). It seems to be a heavily-reddened 
OB supergiant (Kerkwijk et al. 1995). Its physical association with the SGR is in little 
doubt since it coincides with the finely-measured VLA position of the plerion's central peak 
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(Vasisht, Frail & Kulkarni 1995). 

Since 1806-20 is the only SGR with a confirmed companion and a radio plerion, 
it is tempting to speculate that the two phenomena are related, and that the plerion 
is accretion-driven. However, SGR 1806-20 has few properties in common with known 
accretion-powered radio nebulae (Vasisht, Frail & Kulkarni 1995). Moreover, the X-ray 
power is much less than the required power in relativistic particles (TD96). Accretion- 
powered outfiows from neutron stars (e.g., jets) can have only a modest outfiow efficiency, 
so the ratio would be the other way around (TD95; TD96). Still, one should keep in mind 
that SGR 1806-20 is much more complex than the other two SGRs. It is possible that 
interaction with, or accretion from, its hot supergiant companion may play some role in 
this system, even if the SGR is a magnetar. More observations are needed. 

If classic GRBs are produced by old halo SGRs, then bursters which remain in binaries, 
like SGR 1806-20, are constrained to not produce too many observable classic GRBs as they 
grow old. If they did, the concentration of GRBs toward the galactic plane could violate 
BATSE limits. In the HBM, this requirement is actually not very restrictive, because 
bursters in binaries have effectively random beaming directions, greatly reducing their 
observable burst rates compared to bursters which become unbound. Furthermore, GRBs 
from old SGRs in binaries might be suppressed by the environmental influences of their 
companions, or because values of Vr that are sufficiently small to avoid binary disruption 
only occur — via the mechanism of §2.4 — in stars with systematically weaker and/or more 
disordered magnetic fields, which penetrate the neutron star core less deeply. Such stars 
would be expected to dissipate their smaller resevoirs of magnetic energy relatively quickly, 
showing less magnetic activity (GRBs) late in their lives. The fact that the detectable 
energy loss from SGR 1806-20 exceeds that of the other SGRs is consistent with this 
conjecture. 
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B. Alfven Wave-Mediated GRBs in Halo Models 

The amplitude 5B of Alfven waves required to to carry a luminosity of 
L^ = 10^^ L41 (AQ^/iii) erg s~^, as appropriate for halo GRBs (cf. eq. |]ll| 



IS 



/■ *^L., V" ,„,^9rV2 {^n-,V'" 



near the neutron star surface (e.g., Blaes et al. 1989, hereafter BBGM). In this equation 
AQ^ = 27r(l — cos(f)b) is the solid angle of GRB emissions, and AQa is the solid angle of 
the neutron star surface through which Alfven waves pass. These Alfven waves presumably 
couple to shear waves in the crust, which have velocity V^ = (fi/ pY^"^, where the shear 
modulus at density p = 10^^ pn gm cm~^ is /i = 8 x 10^^ pj^^'^'' erg cm~^, valid below neutron 
drip, pii < 4.6. This value of p{p) comes from eq. (11) of TD95 and a power-law fit to 
the equilibrium nuclear parameters in the Baym, Pethick & Sutherland (1971) equation of 
state. For a neutron star with (uniform, vertical) magnetic flux density B, the coupling 
between Alfven and seismic waves occurs at the density px where B = [Arc p^p^)]^^"^, or 

p, = 10^^ Blf gm cm-3. (B2) 

Matching the seismic energy flux Fg = pV^e^, where e is the wave strain amplitude deep in 
the crust, to the surface Alfven energy flux F4 = c((5i?)^/47r implies a strain amplitude at 
the coupling depth of 



Since this ignores wave reflection, it should be considered a lower bound on e. This 
exceeds the dynamical yield (fracturing) strain for B ~ 10^^ G. Thus waves with sufficient 
amplitude to mediate a halo GRB would crum,hle a radiopulsar's crust and thereby dissipate 
energy, probably producing an excess of X-rays. Only for fields that satisfy 

5 > 2 X 10^^ L°f '-^) ^) G, (B4) 



Vio-V \An, 
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is the vibration amplitude small enough to avoid pulverizing the crust. This condition is 
numerically similar to the magnetic energy requirement of equation (jl^), but conceptually 
distinct. 



Note that when equation ( P4|) is satisfied, there are actually four distinct physical 
domains for the propagation of waves in the crust: (1) p > pdrip ~ 4.6 x 10^^ gm cm~^. 
At this depth, the shear waves have amplitude that varies with density as e oc p°'^^, 
based on a power-law fit to the Negele & Vautherin (1973) equilibrium composition. 
(2) px < p < Pdrip, where p^ is given by eqn. ( p2D . Shear waves with e oc p^°-^. (3) 
Prei < P < Pxi where prei = (-B^/47rc^) = 9 X lO^-B^g gm cm~^. Non-relativistic Alfven waves 
with e oc p°'^^. (4) p < prei. Relativistic Alfven waves with e = 6B/B = constant. Thus 
e(p^) as given by eqn. ( P3|) is the maximum strain in the crust above the neutron drip 
layer J^ 

Linear Alfven waves in a neutron star magnetosphere are accompanied by oscillating 
currents, 6j = (c/47r)(V x (5B) — (l/47r)(95E/9t). If the ambient plasma density is so low 
that the drift velocities needed to sustain 6j exceed c, then the wave is charge-starved. This 
occurs at plasma densities less than 

n, ~ ^4 — - = 1 X 10^^ 4f z/3 cm-^ (B5) 

2ec 

where u = 10^ z/3 Hz is the wave frequency. [Note that BBGM predict a frequency spectrum 
peaked at z/ ~ 10^ -E41 Hz for bursts triggered by crust fractures. This comes from the 
condition that z/~^ be comparable to the shear wave crossing-time of the fracture.] BBGM 



^^BBGM did a careful study of the coupling of seismic and Alfven waves, including wave 
refiection effects. However, their numerical results for B = 10^^ G are in the regime in which 
zone (3) effectively does not exist (i.e., px occurs at a depth that is less than one wavelength 
for Alfven waves in the relevant frequency range). 
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have speculated that charge-starved Alfven waves accelerate charges until they are limited 
by radiation losses, thus producing GRBs as the waves undergo nonlinear damping. If this 
is true, then GRB emission from ordinary radiopulsars might occur only after they cross 
the "death hue" at age ~ 2 x 10"^ B^2'^^ Y^s (Chen & Ruderman 1993), after which time 
the ambient plasma density in the magnetosphere is less than n^:. This could produce a 
delayed turn-on of GRBs from pulsars, as in the galactic halo model of Li & Dermer (1992). 
However, the energy requirements for bursters (eq. [|I^ and eq. [0]), the crust-crumbling 
limit (eq. ||B4|| ) and the over-concentration of old pulsars to the galactic disk are all potential 
difficulties with this scenario (LD92). Magnetars, which may avoid these drawbacks, cross 
the death line at t/5 ~ 2 x 10^ Bil yrs, comparable to the length of the SGR activity epoch. 
Subsequently, there exists a co-rotation charge density in the magnetosphere of magnitude 
(Goldreich & Juhan 1969) 

^ ^ 1 X 10^2 5I/2 (±y'^' cm-^ (B6) 



-ProteC \t[), 

where we have assumed a value of Prot{t) appropriate for spindown driven by a constant 
dipole field. Comparison with eq. ( pSJ) shows that, just after crossing the death line, Alfven 
waves carrying luminosities characteristic of halo GRBs are not charge-starved near a 
magnetar's surface. The waves become charge-starved only after a time period 

t, = 4 X 10^ Lj/ z/3-2 yrs, (B7) 

independent of B in this high-i? domain. Interestingly, this is close to the turn-on time 
required by LD92. The connection of charge-starved Alfven waves to GRBs is uncertain, 
and merits further study. 
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FIGURE CAPTIONS 

Fig. 1. — Schematic diagram of a beamed gamma-ray burster escaping from the Galaxy. 
The Galaxy is depicted in cross-section by a bulged disk, oriented horizontally at the bottom 
of the figure. The burster is born at point A in the galactic disk, with a recoil velocity Vr 
in a random direction. Gamma ray emissions occur within indicated cones of angular radius 
06 about the star's magnetic axis /i, which we assume is nearly aligned with Vr. As the 
star sails into the halo, reaching point B and beyond, the line of sight from Earth to the 
star gets increasingly aligned with Vr. Thus more and more escaping bursters become 
visible at larger distances r from the galactic center. The fraction of all bursters which are 
detectable increases with distance as the transverse area of the beaming cone in the galaxy, 
or ~ r^. This tends to counter-balance the free-streaming density trend, n oc r~^ within 
a "core radius" Re ~ Ro/(/)b, where Ro ~ 10 kpc is a typical galactic disk dimension. At 
distances larger than R^., all bursters are detectable at Earth, and the n oc r~^ density trend 
prevails. Because only a tiny fraction 1 — cos06 ~ 0.06 (0fe/2O°)^ of all bursters are visible 
when they are born in the galactic disk, the observable dipole anisotropy of GRB positions, 
(cos0), is greatly reduced compared to what it would be if bursts were emitted isotropically. 
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Fig. 2. — The 2-D cross section of the "zone of avoidance" (ZOA) for beaming angles 
0{, = 10°, 20° and 40°. The thick horizontal line represents the galactic plane where the 
galactic center (GC) and Earth are located. The circles are loci within which bursting stars 
born at GC and moving in straight lines are not detectable at Earth. Everywhere within 
these ZOA circles, the angles between the star-to-Earth and star-to-GC directions are larger 
than 0fe. A star becomes detectable either when it travels beyond the circle, or when it is 
inside the lens-shaped overlap region of circles directly between GC and Earth. The radius 
of the circles are Dsun/siii06, where -Dgun = 8-5 kpc. 

Fig. 3. — Parmeters describing the spatial distribution of detectable bursters in the HBM 
are plotted as a function of distance from Earth, for several different beaming angles 0^: 
90°, or unbeamed (thick solid line); 40° (thin solid line); 20° (dotted); 10° (dashed) and 5° 
(alternating short and long-dashed). The thick long-dashed horizontal lines in each subplot 
are the expected values if the distribution is isotropic at all distances. Filled circles, placed 
at D = 180 kpc in each subplot, are current BATSE values for 1122 bursts; the error bars 
show 3-0" statistical uncertainties. Models with 0f, > 40° or 05 < 10° are evidently ruled out 
by current observational limits, but a significant range of BATSE sampling depths Db ~ 100 
kpc are allowed if 0f, ~ 20°. 

Fig. 4. — Comparison of the intensity distribution of the HBM with combined PVO and 
BATSE data. Note the homogeneous, Euclidean "core" with slope —3/2 at the bright end, 
which turns over to a significantly milder slope at the faint end. A beaming angle 0{, = 20° is 
used in the HBM plotted here, and the (standard-candle) intrinsic burst luminosity is varied 
to make the fitting, which has a x^ ~ 6.9 with 12 degrees of freedom. The implied peak 
luminosity, including the beaming reduction factor, is ~ 1.5 x 10^° ergs s^^. 
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Fig. 5. — Comparison of the intensity distribution of HBM with BATSE only. Again, 
(pb = 20° but the derived peak luminosity from this fitting is slightly smaller than that in 
Figure ^. 

Fig. 6. — Cumulative plots of the angular (an)isotropy measures versus peak photon flux 
P256) for bursts in the 3B catalog (814 in total). The solid curve is the HBM prediction 
for Vr = 1000 km s^^ and (pb = 20°; data points are from BATSE; the dashed line shows 
expected values if the bursts are distributed isotropically on the sky. The intrinsic burst 
peak luminosity found by fitting the log A/" — log P distribution was used to fix the the 
HBM curves' horizontal positions. Error bars on the BATSE points show l-a statistical 
uncertainties appropriate for comparing the data with the HBM. Both models (HBM and 
isotropy) are corrected for BATSE's imperfect sky coverage. Isotropy gives a better fit to 
the data than does the HBM. However, the data are consistent at the ^ 2-cr level or better 
with both models, at all observed fiux levels. 

Fig. 7. — Gamma-ray intensity, I{0), is plotted as a function of angle 6 away from the beam 
axis, for the relativistic beam model of Mao & Yi (1994). In this model, luminous material 
is expanding uniformly into a cone of half opening-angle A6 with Lorentz factor 7, and 
emitting photons isotropically in the material rest frame. In this plot we keep A^ = 20°, 
while varying 7 as shown in the figure legend. Also shown is the (pb = 20° uniform beam 
case (solid line), which is equivalent to A6 = 20° and 7 = cxd. 

Fig. 8. — The amounts by which observable GRB angular statistics in galactic halo models 
with variable beaming patterns I (6) differ iiom their values in the uniform beam ( "unibeam" ) 
HBM, (pb = 20°, are plotted as a function of peak photon fiux. The beam intensity 
models considered here have A9 = 20° and 7 = 3,5, 10 as indicated in the figure legend; 
corresponding I{6) functions are shown in Figure 7. The curves are not completely smooth 
because of the finite number of stars used in the Monte Carlo simulations. 
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Fig. 9. — Similar to Figure || but including the effects of luminosity functions. Filled dots are 
from BATSE; 3-cr error bars are shown at the faint end. The thick solid line is for the HBM 
assuming mono-luminosity GRBs (standard candles). The broken lines are HBM results for 
two alternative luminosity functions discussed in the text: (piL) oc L^^ for Lmin < L < Lmax 
with a dynamic range Lmax/-^min = 20 [dotted line]; and (piL) oc L^ with -Lmax/-^min = 4 
[dashed line]. It is clear that a GRB luminosity function with a positive (negative) slope 
tends to improve (worsen) the HBM fit to BATSE angular data. 
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